Three-state mean-field Potts model with first- and second-order phase transitions 
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We analyze a three-state Potts model built over a ring, with coupling Jo, and the fully connected 
graph, with coupling J. This model is an effective mean-field and can be solved exactly by using 
transfer-matrix method and Cardano formula. When J and Jo are both positive, the model has 
a first-order phase transition which turns out to be a smooth modification of the known phase 
transition of the traditional mean- field Potts model (J > and Jo = 0), despite the connected 
correlation functions are now non zero. However, when J is positive and Jo negative, besides the 
first-order transition, there appears also a hidden (non stable) continuous transition. When J is 
negative the model does not own a phase transition but, interestingly, the dynamics induced by the 
mean-field equations leads to stable orbits of period 2 with a second-order phase transition and with 
the classical critical exponent (3 — 1/2, like in the Ising model. 

PACS numbers: 05.50.-|-q, 64.60.aq, 64.70.-p 



I. INTRODUCTION 

The mean-field concept is one the fundamental 
paradigm in all theoretical physics and its interdisci- 
plinary branches and applications. The mean- field idea 
consists in replacing the interactions acting on a par- 
ticle with an effective external field to be determined 
self-consistently. This paradigm allows one, not only 
to analytically face, in a first approximation, any given 
model, but it also provides a powerful understanding of 
the physics of the model which, even though approxi- 
mate, is often pedagogically deeper than the understand- 
ing one would get from a possible exact solution (if any). 
In particular, it would be harder to understand the con- 
cept of the collective behavior and the phase transitions 
of a system without a suitable mean-field theory. 

At the mathematical base of the mean-field theory 
there are models which are exactly solvable by a mean- 
field technique: the mean-field models. These models 
represent the limit cases of more realistic models in which 
one or more parameters are typically send to or to oo so 
that the mean-field approximation becomes exact. Tra- 
ditionally, the concept of the mean-field models is asso- 
ciated with the absence of correlations in the thermody- 
namic limit. In [fl] (see also [3] and [i]) we have shown 
that this condition is only a sufficient condition for the 
system to be mean-field, but in general it is not necessary. 
There are in fact infinite models having non zero correla- 
tions but still they are mean-field. For example, if Hq is 
an arbitrary Hamiltonian, the model H = Hq+AH, with 
AH a generic fully connected interaction, is mean-field, 
in the sense that we can exactly replace the interactions 
acting on a particle with an effective external field to be 
determined self-consistently. However, now, the presence 
of the term Hq gives rise to non zero correlations when- 
ever Hq has short-range interactions fl5'|. 

Of course, unlike the traditional mean-field models 
(where Hq — and there are not short-range correla- 
tions), the arbitrariness of Hq lets it open now a very 



richer scenario of phase transitions. In particular, it can 
be shown that, when Hq has negative interactions, inver- 
sion transition phenomena and first-order phase transi- 
tions may set in Q . More in general, the phase transition 
scenario associated to the term AH or to H can change 
drastically when Hq has negative couplings. 

In recent years, a renewed attention towards mod- 
els having both short- and long-range interactions, has 
been drawn due to the the importance of small-world 
networks 0, where a finite-dimensional and an infinite- 
dimensional character are both present in the network 
structure. As expected, such models turn out to be 
mean-field, at least for what concerns their critical be- 
havior. However, rather than a theorem, except for the 
Ising case near the critical point 0, and a few ex- 
amples in one dimension @, 0], this turns out to be an 
empirical fact. An exact analytically treatment, even not 
rigorous and confined to relatively simple models, is still 
quite far from being reached out when short-range corre- 
lations are present, as happens in a generic small-world 
network. On the other hand, the models introduced in 
jlj can be seen as ideal small-world networks in which 
the random connectivity of the graph goes to the system 
size N and the coupling J associated to the long-range 
interactions is replaced by J/N. Clearly, without a se- 
rious understanding of the more basic models presented 
in the analytical study of the small- world networks 
and its generalizations (including the scale-free case Q, 
which for the Ising case has been analyzed in [IS]), will 
remain impossible. 

In this spirit, in the present paper we analyze a sim- 
ple and yet rich model: a mean-field Potts model in 
which Hq is a one-dimensional three-state Potts model 
and AH is the traditional three-state mean-field term, 
i.e., the ordinary fully-connected interaction. The mean- 
field equations in this case are sufficiently simple to be ex- 
actly and explicitly solved via the transfer matrix method 
and the Cardano formula for cubic equations. As ex- 
pected, the presence of a non zero positive coupling. 
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Jo > 0, in Hq, alters only smoothly the phase diagram 
of the system (but not its correlations which unlike the 
traditional mean- field term will be now non zero), charac- 
terized by a first-order phase transition. However, when 
Hq has a negative coupling, Jq < 0, besides the first- 
order transition there emerges a new transition which, 
quite surprisingly, is continuous. This transition is not 
stable (the corresponding free energy being not a local 
minimum but a saddle point), so that it does not have a 
physical ground, however it remains a non trivial solution 
of the mean-field equations with a corresponding ground 
state energy (the solution exists at all temperatures) . 

Finally, we analyze the situation in which the term 
AH has a negative coupling J < 0. In this case the 
mean-field equations have only the trivial paramagnetic 
solution and there is no phase transition. However, the 
mean-field equations lead to a dynamics which turns out 
to be stable with an attractor of period 2. Quite inter- 
estingly, this dynamics gives rise only to a second-order 
phase transition with classical critical exponent 1/2, like 
in the Ising model. 



II. GENERAL MEAN-FIELD POTTS MODELS 



in the thermodynamic limit, the order parameter for H, 
Xa — {S{ai,a)), satisfies the system 

Xa = XQ.a{l3Jxi,...,l3JXq) , a = l,...,q (4) 

and the free energy / is given by 

f3f = J2^+Pfoif3Jxi,....PJx,). (5) 

Finally, from Eqs. Q and ([5]) written for an arbi- 
trary external field, it is easy to check that for the 
connected correlation function C-j = '^si^<^i,'^^<^,<^]) ^ 
Y.ai^'^^,'')i^'^:'^j) it holds 

= Col (/^-^^i ' • • • ' /^-^^^ ) + O (^) > (6) 

where Cq.-j is the connected correlation functions for Hq 
in the thermodynamic limit, and the last term in the 
rhs is a finite size effect that can also be calculated 
However, for the sake of simplicity in this paper we will 
not calculate the correlation functions. 



In the spirit of [1], we introduce now a model built 
by using both finite-dimensional and infinite dimensional 
Hamiltonian terms. A generic mean-field Potts model, 
i.e., a model where each variable cr can take q values, 
a = I, . . . , q, can be defined through the following Hamil- 
tonian 



J 



(1) 



where 5((T, a') is the Kronecker delta function and is a 
generic q-states Potts Hamiltonian with no external field. 
Let us rewrite H as (up to terms negligible for N -> oo) 



H = Ho{{a,})~^Yl 



N 



^5((T,;,cr) 



(2) 



As done in [li, fromEq. ^ we see that, by introducing q 
independent Gaussian variables Xa, we can evaluate the 
partition function, Z — J^im} ^^Pi~P^iWi}))j 



Z cx 



-+l3fo(l3Jxi,...,l3Jx,) 



(3) 



where /o(/3/ii, ■ • ■ , (3hq) is the free energy density of the 
Potts model governed by Hq at the temperature 1//3 
and in the presence of a g-component external field 

def 

h — {hi, . . . ,hq). By using the saddle point method, 
from Eq. ^ we find that, if XQ.rj{/3hi, . . . f3hq) is the 
order parameter for Hq as a function of a g-component 
external field [16], {6{ai,<7))Q = XQ-a-iPhi, . . ./Shq), then. 



III. THE TRADITIONAL MEAN-FIELD POTTS 
MODEL 

Before facing the analysis of our model, we want to 
briefly recall the traditional mean-field Potts model de- 
fined as in Eq. ([T]) with Hq = 0. 



A. The pure model 

The use of Eqs. (01)- ([S]) in this case may seem not 
necessary but it is instructive. To apply Eqs. dH)-® 
to the present case, we need to solve the corresponding 
pure model, which is a null Potts model in the presence 
of a uniform external field, h, i.e., we have to calculate 
the following trivial partition function, ZQ{h), which dif- 
fers from Z for the absence of the fully-connected (long- 
range) interaction: 



(7) 



where iV, 



^Yj^5{<7,fJi). We have 
XQ-a {I3hi,...,l3hq) 



^/o(/3/ii,...,/3^) = -logK], 



(8) 



(9) 
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FIG. 1: (Color online) Magnetizations and free energy for the 
case J > and Jo — ^ (equivalent to the traditional mean- 
field Potts model). 



B. The mean-field model 

By plugging Eqs. in Eqs. we get im- 

mediately the following system of equations and the free 
energy density: 



-, CT=1,. 



/3/ = -logK]e''^- +^ 



(10) 



(11) 



To apply Eqs. ((IJ-dS]) to our case we need to solve 
the corresponding pure model, which is a ID three-state 
Potts model in the presence of a three-component uni- 
form external field, h, i.e., we have to calculate the fol- 
lowing partition function 

Zoih)= J2 e^'^°^"i^(^-^-+i)-''^>'^'. (14) 



For any finite N, we can express as ("transfer matrix 
method" ) 



Zoih)=TTT'\ 



(15) 



where T is the 3x3 matrix whose elements, T{(t, a'), for 
any a, a' £ {1, 2, 3}, are defined as 



r(cr, cr') — exp 



pjoS{a,a') + ^(3h, + ^/3ha 



(16) 



For the thermodynamic limit it will be enough to evaluate 
the eigenvalues of T, Ai , A2 , A3 , the free energy density 
of the pure model, /o, being given by 

- l3U/3h) = lim = log(A„,,), (17) 



where A^av is such that I An 



= max{|Ai|, IA2I, IA3I}. 
From (ITBl) we see that the eigenvalues equation reads 

A^ + asA^ + aiA + flo = 0, (18) 

where 

ao = (3e^^° - e^^-'° ~ 2) (e^'^^ + e^''- + e"^^) , (19) 



Eqs. (fTO|) -(fTT |) give rise to a well known phase transition 
scenario 12]: a second-order mean-field Ising phase tran- 
sition sets up only for q = 2, while for any 5 > 3 there 
is a first-order phase transition at the critical value (see 
Fig. [J): 

/3,J=fcillog(g-l). (12) 



N 



(13) 



where we have assumed periodic boundary conditions 
aN+i = o"!, and from now on it is understood that each 
Potts variable a can take the values 1,2, and 3. 



a2 = -e2/3J0 (e^'H+e^''2+g,3h3) 



(21) 



Eq. is cubic in A so that we can solve it explicitly 
using the Cardano formula which gives the three roots: 



IV. THE THREE-STATE ID CASE 

We now specialize the above general result to the case 
in which Hq represents a one-dimensional three-state where 
Potts Hamiltonian: 



Ai = -^a2 + iS + T), 

A2 = -^«2-i(^ + T) + ^(5-T), 

A3 = -^a2-i(^ + T)-^(5-T), 



S 



T 



D = Q^ + R^, Q 



R- 



3ai — 02^ 



9 



R 



9aia2 — 27ao — 2a2^ 
54 ■ 



(22) 
(23) 
(24) 

(25) 

(26) 
(27) 
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Notice that Ai, A2, and A3 are functions of the vector- 
field h. Once the eigenvalues have been calculated, the 
magnetizations, xo-_a-, of the pure model in the thermo- 
dynamic limit can be calculated from Eq. (jl7p as 



1 dX„ 



dfiha 



(28) 



From Eqs. p ^ - (P7|) we see that, in order to evaluate 
Eq. (pS)) . we need to take into account the following 
derivatives, with {a, a', a"} = {1,2,3} 



dap 



(29) 



9a 1 
Who 



dR 



1 dD 



dT 



5 («--*)"'( 



dR 



1 dD 



dD 



ofiK 



+ 2R 



dR 

dWa 



dQ 

dPK 



) dai 



2a2 



aa2 



(31) 



, (32) 



, (33) 



(34) 



(35) 



to Eqs. (2])- ([5]), the role that a generic external field had 
on the pure model Hq, has been now replaced by the ef- 
fective magnetizations Jxc to be found self-consistently 
by Eqs. ([57]) . We could easily consider the presence of an 
additional external field h by simply performing the effec- 
tive substitutions h„ — Jx^+h^ in Eqs. (jl7p - (|36p . which 
does not change the structure of the self-consistent Eqs. 
([57)) , but its numerical detailed analysis goes beyond the 
aim of the present work. 



V. PHYSICAL AND NUMERICAL ANALYSIS 
OF THE SELF-CONSISTENT EQUATIONS 

In this Section we analyze numerically Eqs. (1371) and 
((38)) and provide the corresponding physical explanation. 
It turns out that Amax coincides always with Ai. We 
find it convenient to distinguish the cases Jo > and 
Jo < both for J > 0. Later on we will consider also 
the case J < 0. As is known, the Potts model, like any 
other short-range model, in one dimension does not own 
a spontaneous symmetry breaking. However, from Eqs. 
(|37| and (|38| we see that, for any positive value of J, the 
model governed by H turns out to be a mean-field model 
so that a phase transition is always expected. Therefore, 
in our numerical experiments, it will be enough to keep 
the value of the long-range coupling fixed at J = 1 and 
to obverse what happens by changing the short-range 
coupling Jo. 

In general, the trivial and symmetric solution xi = 
X2 = X3 = 1/3 is always dynamically stable even below 
the critical temperature Tc within a finite range of tem- 
peratures [Ti,Tc], though in general is not leading, and 
below Ti becomes unstable. 



The case J > 0, Jo > 



dR 



«2 + 9^ai-27^-6a 



dag 



dl3h„ 



,2 da2 
'2 9/3 



a2. 



54 



The mean-field model 



-(36) 



By performing the effective substitutions ha — > Jxa in 
Eqs. HDl-dM]), Eqs. ^-§^1 take the form 



dXm.. {{/3h'^}) 



Amax({/3Ja;;}) d/Sha 



\{h'^=Jx'j, (37) 



(3Jxl 



log(A,„ax({/3Jx;})), 



(38) 



where we have written the explicit dependence on the 
arguments {/SJx'^}. 

Note that in the present work we are not going to con- 
sider an additional external field (see Eq. (ID)): according 



When we set Jq — > 0, our model coincides with the 
traditional mean-field model governed by Eqs. ()10p and 
dTTI) . Note however that Eqs. ()18p - (|5S)) are singular for 
Jo = 0. Therefore, to analyze the case Jq = we can 
simply choose a very small value of Jo. From Eq. ()12p 
with g = 3, we see that for J = 1 a first-order phase 
transition develops at the critical point Tc — 0.3607, see 
Fig. ([Ij. The phase transition is triggered by a broken 
symmetry mechanism according to which one of the three 
components (1, 2, 3) becomes favored in spite of the other 
two that remain equal to each other so that, for T — )■ 0, 
two components go to and the favored one reaches the 
value 1. As expected, when Jq > we observe a smooth 
modification of such a scenario, as reported in Figs. ([2])- 
([5]). However, while for relatively small values of Jo we 
see an enhancing of the phase transition, with an increase 
of the corresponding Tc( Jo) as Jo increases, it turns out 
that, for large enough values of Jq, with 0(Jo) — 0{J), 
J'c(Jo) ends to be a monotone increasing function of Jo. 
We note however that, for Jq » J, Tc{Jo) reaches an 
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FIG. 2: (Color online) Magnetizations and free energy for a 
case J > and Jo > 0. 



FIG. 3: (Color online) As in Fig. (2). For visual convenience 
the plots for / are shifted by an irrelevant constant. 



asymptotic value not far from Tc{0). This scenario can 
be partially physically explained by the observation that, 
from one hand, for any J > the model is effectively 
mean-field, so that Tc > but, on the other hand, when 
Jo 3> J, the model remains partially dominated by the 
one-dimensional geometry, and in one dimension (i.e., in 
the pure model) there is no phase transition. The conse- 
quence of these two features is a compromise in which T^, 
remains "small" but finite and, for Jq — >■ oo, Tc reaches 
an asymptotic value of the order of Tc{Q). Note however 
that in the Ising case such an argument does not work: 
when Jo — ?> oo also Tc —J' oo Q. The difference is due to 
the fact that, in the Ising case, the possible transition is 
second-order and regulated by an equation involving xo, 
the susceptibility of the pure model, which has always the 
property xo — ^ oo for Jo — )■ oo , while such an equation 
does not apply for a first-order phase transition, whose 
Tc is not regulated by xo- 
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FIG. 4: (Color online) As in Fig. (2). For visual convenience 
the plots for / are shifted by an irrelevant constant. 



2. The case J > 0, Jo < 

When Jo < the model is still mean field, so that a 
first-order phase transition similar to the case Jq > is 
still present, but with a corresponding lower value for 
Tc(Jo). Now, however, due to the fact that Jo < 0, 
we must to take into account that two spins that are 
consecutive along the ID chain, at low enough tempera- 
tures, cannot have the same value, so that, among the 
three components (1,2,3), the favored one(s), if any, 
along the ID chain must be alternated with another one 
(others). There are two ways to realize this alternation. 
If for example we look for situations in which, at low 
enough temperatures, the component 1 is favored, along 
the ID chain we can look for configurations of the kind 
(1, 2, 1, 3, 1, 2, 1,3,.. .). But we can also look for situa- 
tions in which, for example, both the components 1 and 



2 are favored and, at low enough temperatures, the con- 
figurations are of the kind (1, 2, 1, 2, 1, 2, 1,2,...). Again, 
even if in one dimension (i.e., in the pure model) there 
is no phase transition, the presence of a finite positive 
J guarantees the existence of such phase transitions, as 
confirmed by Figs. (|6]) and 0. Notice, for both the 
situations, with respect to the case Jo > 0, the neces- 
sary modification of the values (xi, X2, X3) for T — ^ due 
to the alternations: now the asymptotic values are ei- 
ther (1/2,1/4,1/4) or (1/2,1/2,0) for the above former 
and latter case, respectively. Quite interestingly, while 
the latter phase-transition mechanism corresponds to a 
first-order phase transition, which turns out to be, in 
shape, quite similar to the phase transition that occurred 
for Jo > 0, the former phase-transition mechanism cor- 
responds to a second-order phase transition. However, 
only the state coming from the first-order transition is 
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FIG. 5: (Color online). As in Fig. (2). For visual convenience 
the plots for / are shifted by an irrelevant constant. 



stable, while the other turns out to be unstable, as seen 
from the fact that the initial conditions giving rise to 
the second-order phase transition live in a subspace of 
(xi,a;2,a;3) of the kind xi — X2 which has zero volume 
in 3 dimensions. A basin of attraction of zero volume 
corresponds to a unstable state. In fact^a control of the 
Hessian of the Landau free energy p7| gives, for the 
solution corresponding to the second-order phase transi- 
tion, always one negative eigenvalue. However, the pres- 
ence a second-order phase transition in a non disordered 
three-state mean-field Potts model, even if hidden in a 
subspace, is a quite non obvious and interesting fact: in 
the subspace xi — X2, as the temperature is decreased, 
the system is forced to favorite X3 not by a jump, but 
continuously. 



A. The case J < 

When J < 0, the approach we have used via Eq. ^ 
is not valid since the Gaussian integral diverges. Yet, 
the saddle point Eqs. (|37|) are still exact as derived from 
the general theorem presented in [l[ (while the free en- 
ergy has a different form with respect to Eq. ([35])). When 
J < 0, Eqs. ([57| have only the trivial symmetric solution 
and no phase transition sets in. However, when analyzed 
dynamically under the iteration Xn+i — F{xn), where 
F{xn) comes from the rhs of Eq. (1571) . the saddle point 
equations develop a phase transition. For J < 0, a part 
from the trivial solution (the paramagnetic/symmetric 
one), we observe in fact that, even if Eqs. p7p do not 
have solution, and the corresponding iteration do not 
converge to any stable point-like solution, Eqs. ([571) still 
own stable trajectories of period 2. In the language of 
dynamical systems we say that the dynamics induced by 
Eqs. ([57|) has an attractor. In the present work, for the 
sake of simplicity, we limit ourselves to the analysis of 



FIG. 6: (Color online). Here there are two unstable states: 
the Unstable A state is obtained by starting with any initial 
condition of the kind X2 = X3 < xi, while the Unstable S state 
is obtained only with the initial condition xi = X2 = X3. Fi- 
nally the leading state is obtained with any initial conditions 
of the kind xi ^ X2 > X3. 
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FIG. 7: (Color online). As in Fig. (6). This case represents 
effectively the limit Jo — >■ — cxj. 



the case Jo — of Eqs. ([57)) . In the Ising case also, the 
dynamics led by the saddle point Eqs. has an attractor of 
period 2 due to the simple fact that mo(— /i) = — mo(ft.), 
where mo(ft.) is the magnetization of the pure model in 
the presence of the external field h. Less trivial is the 
Potts case: it is not analitically obvious what is the re- 
lation of the rhs of Eqs. (|37p under the change J — — J. 
When seen within the iterative dynamics, see Fig. ([8|), 
the model ([T]) presents a phase-transition with a Tc which 
is of the same order of magnitude of the case J > but, 
quite interestingly, the phase transition is second-order, 
with an Ising-like shape having the classical mean-field 
critical exponent for each component /3 = 0.5 (not to be 
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FIG. 8: (Color online) Stable magnetizations of the period- 
2 attractor for the case J = — 1. In the inset enlargement 
around the transition point where, for T < Tc, we plot also 
the matching functions y{T) = 1/3 + 0.65(rc - and 
z{T) = 1/3 + 0.35(rc - T)^/^ with = 1/3. 



confused with the inverse of the critical temperature), as 
can be seen from the excellent agreement with the match- 
ing functions y{T) and z{T) plotted in the enlargement 
of Fig. ([5]). A similar behavior is expected for Jg ^ 0. 

The iterative dynamics Xn+i = F{xn), or any other 
variations of it, has not, a priori, any direct relation with 
the real model. However, the Glauber dynamics 13] of 
the model ([T]), at least in its discretized version, repro- 
duce a similar behavior since the stationary equations 
of the Glauber dynamics have the same structure of the 
equations found in our Gibbs-Boltzmann thermodynamic 
context. Alternatively, one can expect a similar behavior 
in a continuous time dynamics when synchronous spin- 
flip updates are considered jl4{ . 



Jo, and the fully connected graph, with coupling J. This 
model is a non trivial exactly solvable mean-field where 
new phenomena emerge as a consequence of the presence 
of a finite-dimensional geometry, where correlations are 
not negligible. For J > we have found, in particular, 
the existence of a finite critical temperature in the limit 
Jo — > oo, which is not seen in the Ising case, and the 
existence of a hidden continuous phase transition for the 
case Jo < 0. Of course, for Jo < the ground state of the 
system has a totally different symmetry with respect to 
the case Jo > (compare the asymptotic values toward 
T = in Figs. (1-7)), like in the pure model (J = 0), but 
with the important novelty that now (J > 0) the sys- 
tem is an effective mean-field model with a finite critical 
temperature. 

In the case J < there is no phase transition, yet 
the dynamics induced by the saddle point equations has 
an attractor of period 2 to which the system converges. 
Quite interestingly, within this dynamics, the system un- 
dergoes a second-order phase transition with the univer- 
sal critical exponent f3 — 1/2, like in the Ising model. We 
think that this phenomena can be similarly found in a 
continuous synchronous time dynamics like, the Glauber 
dynamics with multi-spin fiip updates, where the station- 
ary equations have the same structure of the equations 
found in our Gibbs-Boltzmann thermodynamic context. 
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